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Abstract 

Background: Mitochondrial 2-oxoglutarate (a-ketoglutarate) dehydrogenase complex (OGDHC), a key regulatory 
point of tricarboxylic acid (TCA) cycle, plays vital roles in multiple pathways of energy metabolism and biosynthesis. 
The catalytic mechanism and allosteric regulation of this large enzyme complex are not fully understood. Here 
computer simulation is used to test possible catalytic mechanisms and mechanisms of allosteric regulation of the 
enzyme by nucleotides (ATP, ADP), pH, and metal ion cofactors (Ca^^ and Mg^^). 

Results: A model was developed based on an ordered ter-ter enzyme kinetic mechanism combined with con- 
formational changes that involve rotation of one lipoic acid between three catalytic sites inside the enzyme 
complex. The model was parameterized using a large number of kinetic data sets on the activity of OGDHC, and 
validated by comparison of model predictions to independent data. 

Conclusions: The developed model suggests a hybrid rapid-equilibrium ping-pong random mechanism for the 
kinetics of OGDHC, consistent with previously reported mechanisms, and accurately describes the experimentally 
observed regulatory effects of cofactors on the OGDHC activity. This analysis provides a single consistent 
theoretical explanation for a number of apparently contradictory results on the roles of phosphorylation potential, 
NAD (H) oxidation-reduction state ratio, as well as the regulatory effects of metal ions on ODGHC function. 



Background 

The 2-oxoglutarate (a-ketoglutarate; aKG) dehydrogen- 
ase complex (OGDHC, EC 1.2.4.2, EC 2.3.1.61, and EC 
1.6.4.3) is a multi-enzyme complex which catalyzes the 
chemical reaction: 



aKG^ - -I- CoASH"" -i- NAD" 



+ H2O ^ 



Succinyl - CoA'*" -1- CO^" -1- NADH^- -1- H + 



(1) 



OGDHC is primarily located within the mitochondrial 
matrix and is a key regulatory enzyme complex in the 
TCA cycle, responsible for oxidative decarboxylation of 
2-oxoglutarate, transferring a succinyl group to coen- 
zyme A (CoASH*') and producing reducing equivalents 
(NADH ') for the electron transport system. Regulation 
of OGDHC not only affects the distribution of 2-oxoglu- 
tarate between the TCA cycle and malate-aspartate 
shuttle system, but also has effects on the oxidative 
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deamination of glutamate. OGDHC is a crucial target of 
reactive oxygen species (ROS) and also able to generate 
ROS, which make it distinctly important for bioener- 
getics [1]. The molecular organization of OGDHC is 
similar to that of the pyruvate dehydrogenase complex 
(PDHC) as it belongs to the same heterogeneous family 
of 2-oxo acid dehydrogenase multi-enzyme complexes 
[2]. It consists of multiple copies of three enzyme com- 
ponents: oxoglutarate dehydrogenase (El), dihydro- 
lipoamide succinyltransferase (E2), and dihydro-lipoa- 
mide dehydrogenase (E3). Consecutive actions of these 
enzymes catalyze the oxidation of 2-oxoglutarate and 
reduction ofNAD', which results in the production of 
NADH^' and Succinyl-CoA*' (Figure lA). Allosteric 
interactions associated with the El component are 
known to be the predominant target for controlling of 
OG-DHC activity [3]. 

OGDHC was first purified from the pig heart mito- 
chondria by Sanadi et al. [4] and subsequently studied 
by many researchers to examine its catalytic and regula- 
tory properties within permeabilized, un-coupled, and 
coupled mitochondria from a variety of mammalian 



o 



© 201 1 Qi et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons 
BiolVlGCl Central Attribution License (http://creativecommons.0rg/licenses/by/2.O), which permits unrestricted use, distribution, and reproduction in 
any medium, provided the original work is properly cited. 



Qi et al. BMC Biochemistry 201 1, 12:53 
http://www.biomedcentral.eom/1 471-2091/1 2/53 



Page 2 of 1 5 



(A) 



AKG CO 



NADH 




NAD 



CoA 



succinyl-CoA 



\ J, LD (Lipoyl Domain) 
BD (Binding Domain) 



Lip 
/\ 

s S 



(B) 



A 



A 



SH S-succinyl 



SH 



SH 




fA 



f-A, 




Figure 1 Schematic representation of tlie proposed mechanism of 2-oxoglutarate dehydrogenase complex (OGDHC) (A) It consists of 
three component enzymes: oxoglutarate dehydrogenase (El), dihydro-lipoamide succinyltransferase (E2), and dihydro-lipoamide dehydrogenase 
(E3). The schematic representation here does not describe true stoichiometry of the multiple copies of three enzymes in the complex. The 
binding domain and lipoyl domain of E2 polypeptide are connected to the complex core with flexible links (dotted lines), are used here to 
describe the mechanism in which a single lipoic acid rotates between the three catalytic sites. In the catalytic cycle, the disulfide at the tip of 
the lipoyl can be in oxidized, reduced or semi-reduced lipoate forms, the later one is connected with succinyl residue transferred from 
oxoglutarate. (B) This schematic illustrates the proposed kinetic schemes along with the mechanism of conformational changes. The forward 
reaction is read in the clockwise direction. The complex has three binding sites: one for each composite enzyme. El binds to 2-oxoglutarate or 
corresponding product CO2 (top), E2 binds to COA or corresponding product Succinyl-CoA (bottom-right), E3 binds to NAD or corresponding 
product NADH (bottom-left). It is assumed that, in the process of conformational changes, the rotation of one lipoic acid between three catalytic 
sites leads to transfer succinyl from El to E2 and proton from E2 to E3. 
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tissues [5-11]. A catalytic mechanism for the overall 
reaction of the enzyme complex was also first proposed 
by Sanadi et al. [4] which suggested that the coenzyme, 
NAD', and 2-oxoglutaric acid participate in the reaction 
with the help of the cofactors thiamine pyrophosphate 
(TPP), lipoic acid, and FAD^' [12,13]. Their proposed 
mechanism is a Hexa-Uni-Ping-Pong mechanism in Cle- 
land's terminology [14] where it is assumed that the first 
product (CO2) is released before the second substrate 
(CoASH* ) binds, and the second product (Succinyl- 
CoA ') is released before the third substrate (NAD ) 
binds to the enzyme. Subsequently, Koike et al. [15] 
postulated another mechanism in which, the lipoic acids 
transfer intermediates by rotating between the three cat- 
alytic sites. Furthermore, experimental results of fluores- 
cence resonance energy transfer and dynamic anisotropy 
showed that the lipoic acids in the E2 component 
undergo motion where they rotate between different 
catalytic sites [16-18]. The results of steady-state kinetic 
studies done by Hamada et al. [19] and Smith et al. [10] 
contradict each other, and not all results are compatible 
with the Sanadi mechanism [4]. This issue was 
addressed by Mcminn and Ottaway [20] with kinetic 
studies based on the Fromm method [21]. Mcminn and 
Ottaway [20] explained the observed nonlinearity in the 
reciprocal plots of the results and proposed a phenom- 
enological mechanism with semi-random characteristic. 
A recent study by Aevarsson et al. [22] on the crystal 
structure and architecture of 2-oxo acid dehydrogenase 
multi-enzyme complexes, provides interesting insights 
into the plausible kinetic mechanism of 2-oxo acid 
dehydrogenase family which includes OGDHC. 

It has been consistently shown that the activity of 
OGDHC is controlled by various factors, including the 
variations of the NAD oxidation-reduction state, the 
state of phosphorylation of the nucleotide systems, and 
the ratio of succinyl-CoA to CoA-SH. Regulation by 
reversible phosphorylation has not been demonstrated. 
Experiments in isolated mitochondria of heart, liver, and 
kidney have shown that the OGDHC is regulated by Ca^ 
^ ions with a marked decrease in the apparent /<,„ for 2- 
oxoglutarate in the presence of adenine nucleotides and 
minimal effect of Ca^^ at saturating concentration of 2- 
oxoglutarate [9,23-25]. Moreover, the apparent for 
2-oxoglutarate is lowered by a decrease in the ATP/ 
ADP ratio, which can significantly increase the sensitiv- 
ity of the enzyme to Ca^^ions [23,25]. It has been 
reported that the maximum activity of OGDHC is unaf- 
fected by changes in pH, while the apparent of the 
enzyme for 2-oxoglutarate is greatly altered by changes 
in pH over the range of 6.5 - 7.5 [23]. A number of stu- 
dies have also demonstrated the possible role of Mg^* 
ions in the regulation of OGDHC either by directly 
affecting the activity of the enzyme or by modulating 



the Ca^^ effect on the enzyme. Mg^^ ion has been 
shown to increase [26-28] or to have no effect on the 
activity of OGDHC [29,30]. McCormack and Denton 
[23] studied isolated OGDHC from pig heart mitochon- 
dria and found that there is no effect of EDTA and 1 
mM Mg^"^ on the activity of OGDHC when Ca^"^ con- 
centration was effectively less that 1 |iM. Panov and 
Scarpa [9] concluded that the effects of Mg * and Ca ^ 
ions on the OGDHC activity are additive only at relative 
low concentration of free cations which suggested that 
at high concentrations, each ion may compete each 
other for binding sites. It is also evident that, in the pre- 
sence of low Ca ^ concentration, Mg * ion can strongly 
modify the enzyme's affinities for 2-oxoglutarate and 
NAD' [9]. However, the kinetic mechanisms by which 
these divalent metal ions regulate the properties of 
mitochon-drial OGDHC are not understood. 

Although a number of attempts have been made to 
understand the catalytic mechanisms of OGDHC, both 
experimentally and theoretically, there is no mechanistic 
model that consistently explains the available experi- 
mental data on the kinetics of this enzyme complex and 
adequately describes the regulatory roles of nucleotides 
and other metal ion cofactors (Ca^*, Mg^*, etc.). There- 
fore a mechanistic model of OGDHC is needed to 
understand the orchestrated controlling of OGDHC by 
cofactors inside mitochondria under different physiolo- 
gical conditions. In the present work, a kinetic model of 
OGDHC is introduced to quantitatively understand the 
catalytic properties and regulation of OGDHC, based on 
the observations from a large number of independent 
experimental studies in mammalian tissues. The model 
accurately describes the catalytic properties of this 
enzyme complex observed experimentally, and clarifies 
many contradictory results reported in earlier studies. 

Methods 

In this section, we first present a general kinetic model 
for conformational changes in OGDHC, based on a pre- 
sumed ter-ter enzyme mechanism via substrate channel- 
ing. The model is then used to characterize the kinetics 
of the 2-oxoglutarate dehydrogenase reaction (Equation 
1) and further extended to describe the regulatory roles 
of cofactors, i.e., nucleotides and various metal ions. 
The kinetic parameters of the model are estimated using 
a wide variety of experimental data, available in the 
literature. 

Kinetic scheme for conformational changes in a ter-ter 
enzyme mechanism 

The kinetic equation of the proposed model for 
OGDHC reaction is derived from a ter-ter enzyme me- 
chanism combined with a model of conformational 
changes that represent the rotation of the single lipoic 
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acid between different catalytic sites [15,18]. The deriva- 
tion is inspired by a previously developed model for 
trans-carboxylase [31]. The assumption of the model is 
that the enzyme complex is composed of three sub- 
enzyme (El, E2, and E3), each with one binding site: 
site 1 binds to 2-oxoglutarate {aKG^') or corresponding 
product CO2, site 2 binds to CoASH*' or corresponding 
product Succinyl-CoA*', site 3 binds toNAD' or corre- 
sponding product NADH^' (Figure lA). Furthermore, 
the basic mechanism involves conformational changes, 
where the rotation of one lipoic acid between three cat- 
alytic sites leads to transfer of succinyl from El to E2 
and proton from E2 to E3. In the catalytic cycle, the dis- 
ulfide at the tip of the lipoyl can be in oxidized, reduced 
or semi-reduced li-poate forms, the semi-reduced form 
is bound with succinyl residue transferred from 2-oxo- 
glutarate (Figure IB). 

Each of the six conformational states shown in Figure 
IB can involve any possible binding states associated 
with the enzyme. For example, the first site is either 
empty or bound to 2-oxoglutarate or CO2; the second 
site is either empty or bound to CoASH ' or Succinyl- 
CoA'*'; and the third site is either empty or bound to 
NAD' or NADH^'. Therefore, there are a total of 27 
binding states for each one of the six conformational 
states, which gives rise to 27 x 6 = 162 distinct states in 
the model. Here we denote these 162 states as E]^^, 
where i e {1, 2, 3, 4, 5, 6} represents the index for con- 
formational states, and x e {0, A, P}, y e {0, B, Q} and 
z e {0, C, R} represent the binding states of site 1, site 
2, and site 3. The lower-case e'^^ is used to represent 
the fraction of each state. Therefore, the total fractional 
states can be expressed as 



^ total ~ S ^xyz- 

;ce(0,A,P),}'e{0,B,Q),ze|0,C,R) 



(2) 



We assume rapid equilibrium binding for all 27 bind- 
ing states, implying that the binding processes are much 
faster than the conformational change processes. With 
this assumption. Equation (2) can be written as: 



^total ~ ^free ^ 



(1 + IA]/Kj, + [P]/Kp) 
X {1 + \B]/Kb + [Q]/Kq) 
X {1 + IC]/Kc + [R]/Kr), 



(3) 



where 6^^,^, = ^000 denotes the fraction of free enzyme 

complex that binds to the reactants; Ka' Kb> Ko Kp, Kq 
and Kr are the dissociation constants associated with 
the binding of reactants (A: aKG^', B: CoASH*' and C: 
NAD') and products (P: CO2, Q: Succinyl-CoA*' and R: 
NADH^') to the enzyme complex. Here, we assumed 
that these constants do not depend on conformational 



states of the enzyme complex and the bound reactant at 
one site does not influence the binding reaction at 
another site: all binding interactions are independent of 
one other. These assumptions are necessary to make the 
model tractable, and are validated by comparing the 
model predictions to the available experimental data. 

We define fi as the fractions in conformation state i (i 
e {1, 3, 5}) which can undergo forward conformational 
transformation to state /+1 (lipoate changes in three dif- 
ferent redox forms). Similarly, is the fraction in con- 
formation state /+1 which can undergo conformational 
transformation in the reverse direction (Figure IB). Spe- 
cifically, the binding of 2-oxoglutarate at site 1 is neces- 
sary for transition from conformation state 1 to 2. 
Therefore, we can write 



fi = 

h = 

fs = 



1 + [K\IKa + [P]/Kp 



1 + IA]/Ka 



'-I V 

+ [n/Kp) 



1 + \B]/Kb + [Q]/KqJ' 
\Q]/Kq 
1 + [B]/Kb + [Q]/Kq 
[C]/Kc 
1 + [C]/Kc + [R]/Kr 

m/KR 

1 + [C]/Kc + [R]/Kr 



(4) 



The net turn-over (reaction velocity) for this mechan- 
ism can be expressed 



V 



[Etotal] 



total 1 



total 1 



(5) 



Applying the King and Altman method to the scheme 
shown in Figure IB gives the following expression for 
the net reaction velocity: 



(6) 



+f3f5k2kskiks{ks + fe_6) + f-jf-sikl + fe-2)fe-3fc-4fe-5fe-6 

+/_3/i(fe2 + k-2]k-3k-4kiks + fifsikj + k-2)kik4ksk6 

+flMki + k-4)kik2kiki + /_i/_5(fe4 + *!_4)fe_5fe_6fe_ife_2 

+fif-5{k4 + k-t]k-sk-6k2k3 + f-if-zfsk-ik-3ksk-2(k6 + fe-s) 
+f-ihfsk-ik3ksk4(k(, + fe-s) + /_3/_5/ik_3fe_5fei*!_4(fe2 + fe-2) 
+f-3fif5k-,kikski(k2 + fe-2) + f-if-5hk-ik-skik-i{k-t + ^4) 
+f-5hhk-skikik2(k4 + fe-4) + fihfskikiks(k4k2 + kf,k, + k2kf,) 

V+/_l/_3/-5fe-lfe-3fe-5(fe-2fe-4 + k-ik-s + k-f,k-2) / 



This complex expression can also be obtained using 
our KAPattern package [32], available freely for the deri- 
vation of enzyme rate equations. Substituting the frac- 
tional occupancy distributions as defined in Equation 
(4), we obtain an expression for the reaction velocity in 
terms of the individual rate constants and dissociation 
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constants. The kinetic constants can be also expressed 
in terms of various rate constants. Using the Haldane 
relationship, the velocity equation can be written as: 

„_ V;V,(|A||B|[Cl-[Pl|Q|[R]/fC„,) 

/ V,[A]|B]|C| + K„,cV,.[A][B][R]/Kir \ 

+KmcVr[Al[B| + J<„BVr[A]|C|[Q|/]f„, 
+ifmBVr|A][C]+Jf„pV||A]|Q]|R|/]C,v,/if,, 

+Jf,cJfmBVr[A]|Q]/Jf,, + Jf^Vr|B|[Cl|Pl/% (7) 
+ifmAVr|B||C] + K,aK„cVrmF][R]/KiplKir 
+Jf,V,JfmcVr[B]|R|/K,r + V,[C| [P] |Q]/]fi,/]fip 

+KibK„AV,lC][P]/Kip + V;[P|[Q|[R|/]f„, 

\+i<mBV/|P|[Q]/]f„, + if,„QVf[P]|R|/]f„, + K,npVf[Q][R]/K,J 

where the kinetic constants are defined as: 

Vf = num\jCoefABC, Vr = tium2lCoefpQR, = CoefBcl Coefmc, KmB = CoefAcjCoefj^c. 

Kmc = CoefAB/CoefABC, KmP = CoefQs/CoefpQR, K^Q = CoefpsjCoefpQR, KmR = CoefpQ/CoefpQs, 

Kip = CoefBR/CoefgpR = CoefzclCoefBcp, Kit, = Coefcp/CoefcpQ = CoefAcjCoefAj^Q, 

Kir = CoeffiQ/CaeffiQR = Coe/As/Coe/ABR, Kia = Coefqsl CoefAQR = Coefss/CoefMR, 

Kib = CoefpajCoefspR = CoefcpjCoefscp, Kic = CoefpQ/CoefcpQ = CoefAQ/CoefACQ and 

Kaj = numi/numz. 

Here we used the shorthand notation similar to that of 
Segel [33] where num = ki 1(2 k^ k^ keKpKoKi^ num = k.^ 
L2 k.3 t4 k.5 teKAKsKc, and Coe^tg = KcKpKqKr {ki ^2 ^3 
ke (fci + ^-4)), and so on. Vf and have units of mass per 
unit time per mass of protein. Other kinetic parameters 
associated with the binding of reactants and products have 
the units of concentration (mass per unit volume). 

In the ter-ter biochemical reaction, the fourteen 
unknown kinetic parameters in Equation (7) are related 
to the equilibrium constant Kgq (known) via the follow- 
ing equilibrium relationship: 



num I 



(8) 



/ [P][Q1[R] \ 

V[A][B][C]Aq num2 

KipKiqKir VfKmpKiqKir 

KiaKiijKic VrKiaKmsKic 
VfKipKmQKir ^ VfKipKujKmR 

VrKuiKihKmC VrKmAKihKic 



reducing the total number of independent unknown 
kinetic parameter to thirteen. 

Kinetic model of OGDHC using a ter-ter enzyme 
mechanism 

We apply the above general form of the ter-ter enzyme 
mechanism for the analysis of available expe-rimental 
data on the kinetic of OGDHC to estimate the unknown 
kinetic parameters and to elucidate whether the proposed 
mechanism is able to explain the available kinetic data. 

In the TCA cycle, OGDHC is primarily involved in 
the fifth step for oxidation of acetyl-CoA. The reference 
chemical reaction is given by Equation (1). The corre- 
sponding biochemical reaction is given by 



aKG + CoASH 

+ NAD ^ succinyl - CoA + COj tot + NADH. 



(9) 



Here biochemical reactants, such as aKG, correspond 
to ensemble chemical species, such as aKG H aKG', 
etc. The chemical reaction in Equation (1) is unambigu- 
ously balanced in terms of mass and charge, whereas 
this biochemical reaction is not. In this reaction, the 
reactant CO2, tot represents the sum of aqueous carbon 
dioxide and bicarbonate species {COf~ , HCO3 and 
H2CO3). 

The equilibrium constant for the reference reaction 
can be written as: 



a _/ [Succinyl -CoA*- 1 ICO3 



|[NADtf-||H + | 



[aKC^- ] [CoASH"- 1 |NAD- 



exp 



RT 



(10) 



where ^rG°^^jj^ is the standard Gibbs free energy of 

the reference reaction which is computed using the 
basic thermodynamic data (298.15 K, I = 0.15 M) Usted 
in Li et al. [34] . 

For the development of the kinetic model of OGDHC, 
we assume that the ter-ter enzyme mechanism proposed 
in the previous section along with the conformational 
changes (Figure IB) can explain the observed kinetics of 
OGDHC. Because the kinetic data we used here to esti- 
mate the unknown kinetic parameters were all from the 
initial velocity studies in which only the products 
NADH ' and Succinyl-CoA' were present in the reac- 
tion mediums, the forward flux of OGDHC can be 
reduced from Equation (7) as: 



V;[A][B][C] 



ogdhc 



/[A][B][C]+i^„c[A][B][R]/K„ \ 
+J^mc[A][B] +i^™B[A][C][Q]/K,, 
+KmB [A] [C] + K,cK„,B [A] [QyKiq 
+K„BKt[A][Q][R]/K,^/Kir 

+K,„Amc]+Ki,K„ciBmyKir 

\+K„,BKicKialQ]lR]/Kiq/Kir 



(11) 



where [A], [B], [C], [R] and [Q] represent the concen- 
trations of aKG^', CoASH*', NAD', NADH^', and Succi- 
nyl-CoA*', respectively. This kinetic expression for 
OGDHC reaction contains 8 unknown kinetic para- 
meters. Thus, this expression was used first to estimate 
the 8 unknown kinetic parameters. Using the relation- 
ship Jogdhc/JogdiK ^ exp(AG/i?T) to determine the reverse 
flux [35], we obtain the full flux expression for OGDHC 
as: 
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= - r 

ogdhc ' ogdhc ' o^dhc 
( 



. J+ 

' ogdhc 



I ISuccinyl - CoA^-nCOj^l^ \ 
1 I |NADH^"][H*] 



K°,,,,d,c ([o-KG^- l[CoASH*-][NAD- 



V 



(12) 



The kinetic expressions for the regulatory effects of 
various cofactors on OGDHC (which depend on 10 
additional kinetic parameters) are parameterized in the 
Results section, and are estimated separately (see 
below). Because CO2, tot dependent terms are not 
included in the denominator in Equation (11), CO2, tot 
dependent product inhibition is not accounted for in 
Equation (12). 

Parameter estimation 

The developed kinetic model of OGDHC has 8 adjusta- 
ble parameters for catalytic mechanism and 10 adjusta- 
ble parameters for cofactor regulation (Table 1). 
Parameter values were estimated in a systematic modu- 
lar manner in multiple steps by least-squares fitting of 
the model simulated steady-state flux to the available 
experimental data as detailed in Results section below. 
The FMINCON algorithm in MATLAB (The Math- 
Works, Natick, MA) was used to solve this non-linear 
optimization problem. In addition, sensitivity analysis 
was performed to estimate the sensitivity of the least 
square error to small changes in the optimal parameter 
values. The sensitivity was computed using: 



max(|£*(x,- ± O.lXi) - E*{x,)\) 
0.1E*{xi) ' 



(13) 



where E"' is the least square difference between model 
simulations and experimental data, and is the opti- 
mized value of the f'' parameter. 

Since all kinetic parameters in the model are mea- 
sured relative to species concentration, we performed a 



composition analysis to estimate the concentration of all 
ionic species all experiments analyzed here [36]. 

Results 

Parameterization of basic kinetic mechanism of OGDhHC 

In this section, we present the detailed parameterization 
and validation of the proposed kinetic model based on 
the available experimental data on the kinetics of 
OGDHC, measured in a wide variety of experimental 
conditions. To study the catalytic mechanism of 
OGDHC, McMinn and Ottaway [20] investigated the 
kinetic properties of the OGDHC system, which was 
prepared from fresh pig heart mitochondria. Following 
the method of Fromm [21] in which it was concluded 
from initial velocity studies that the catalytic mechanism 
of OGDHC is not consistent with the Hexa-Uni-Ping- 
Pong mechanism. While their observations suggest a 
random order kinetic mechanism with respect to the 
binding of NAD' and CoASH*' and release of Succinyl- 
CoA*', the binding of 2-oxoglutarate and release of CO2 
is described as a Ping-Pong mechanism. Initial velocity 
kinetics measured by Smith et al. [10] with purified pig 
heart mitochondria OGDHC showed that Succinyl- 
CoA*' and NADH^' were inhibitors, but no inhibitory 
effects were observed with GTP or ATP. Their results 
also show that Succinyl-CoA"*' inhibition was competi- 
tive with CoASH ' and independent of the NAD' oxida- 
tion-reduction state. These data are used here to 
identify the kinetic parameters of our OGDHC model. 

The experimental data in Figure 2 were used to esti- 
mate the values of unknown kinetic parameters that 
govern the basic catalytic mechanism of OGDHC 
(Equation (12)) based on the best fits of the model to 
the data (See Figure caption for details). Measured 
enzyme activity is expressed in [imoles NADH^' 
formed/mg protein/min. We follow a systematic opti- 
mization procedure to estimate each kinetic para- 
meter of the model using appropriate experimental 



Table 1 Kinetic parameter values for 2-oxoglutarate dehydrogenase complex 



Basic kinetic parameters Kinetic parameters for regulatory cofactors 

Parameter: Value Sensitivity Parameter: Value Sensitivity 



Vmay. ^Hiol mg"^ min"''' 






KaCa (MM) 


0.893 


0.293 


KmA (mM) 


0.273 


3.601 


KiATp (mM) 


0.106 


0.943 


KmB (MM) 


6.96 


0.934 


KiADP (mM) 


0.305 


2.073 


Kmc im 


98.6 


2.839 


KaMg (MM) 


19.49 


0.228 


Kia (mM) 


75.9 


1.981 


"Ca 


0.262 


9.648 


Kir (mM) 


2.4 


1.995 


ccajp 


6.694 


1.628 


Ki, (mM) 


0.112 


3.970 


ccadp 


0.173 


4412 


K„ (nM) 


0.218 


5.224 




1.00 


6.594 


KaH (pH) 


6.11 


3.623 




4222 


4.603 



^ We cannot estimate the value of parameter Vf, max in the model without knowing the enzyme activity in the experimental setup. The value of l/max reported in 
the table is equal X x Vf, ^ax- where X is the enzyme activity. The optimization estimates that l/max ~ 2.16 [jmol mg'^ min'^ for Figure 2, \/max ~ T 1.033 unit/mg for 
Figure 3 and Figure 4, Umax ~ 167 nmol mg"^ min"^ for Figure 5. 
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Figure 2 OGDHC activity as a function of different substrate and product concentrations. (A) OGDHC activity as a function of 2- 
oxoglutarate concentration at different NAD and CoASH concentrations in tlie absence of products. Kinetic data were obtained from Figure One 
A of [20]. Concentration of substrates NAD and CoASH in the assay were fixed to: 0.033 mM and 0.005 mM (square); 0.066 mM and 0.01 mM 
(down-triangle); 0.133 mM and 0.02 mM (up-triangle); 0.333 mM and 0.05 mM (circle). (B) OGDHC activity as a function of CoASH concentration 
at different 2-oxog!utarate and NAD concentrations in tine absence of products. Kinetic data were obtained from Figure One B of [20]. 
Concentration of tine substrates 2-oxoglutarate and NAD in the assay were fixed to: 0.025 mM and 0.02 mM (square); 0.05 mM and 0.04 mM 
(down-triangle); 0.1 mM and 0.08 mM (up-triangle); and 0.5 mM and 04 mM (circle). (C) OGDHC activity as a function of NAD at different 2- 
oxoglutarate and CoA concentrations in the absence of products. Kinetic data were obtained from Figure One C of [20]. Concentration of the 
substrates 2-oxoglutarate and CoASH in the assay were fixed to: 0.05 mM and 0.005 mM (square); 0.1 mM and 0.01 mM (down-triangle); 0.2 mM 
and 0.02 mM (up-triangle); and 0.5 mM and 0.05 mM (circle). (D) OGDHC activit/ as a function of NAD at different 2-oxoglutarate and CoASH 
concentrations in the presence of products NAGH. Kinetic data were obtained from Figure Two B of [20] and rescaled for data consistence. 
Concentration of the substrates 2-oxoglutarate and CoASH in the assay were fixed to: 0.5 mM and 0.05 mM. NADH concentration were 0.0 
(circle), 0.01 mM (up-triangle), 0.02 mM (down-triangle) and 0.05 mM (square). (E) OGDHC activity as a function of CoASH concentration in the 
presence of product Succinyl-CoA. Kinetic data were obtained from Figure Five A of [10]. Concentration of the substrates 2-oxoglutarate and 
NAD in the assay were fixed to: 1.0 mM and 336 (iM. Succinyl-CoA concentration were 0.0 (circle), 6.6 pM (up-triangle), and 13.2 |jM (down- 
triangle). (F) OGDHC activity as a function of CoASH concentration in the presence of products Succinyl-CoA and NADH. Kinetic data were 
obtained from Figure Five B of [10]. Concentration of the substrates 2-oxoglutarate and NAD in the assay were fixed to: 1.0 mM and 16 |jM. 
NADH concentration was 42 |jM. Succinyl-CoA concentration were 0.0 (circle), 6.6 ^iM (up-triangle), and 13.2 pM (down-triangle). Solid lines are 
model fitting results from the data points represented by symbols. Experimental data were obtained in 0.08 mM potassium phosphate buffer at 
pH 7.2 and 30°C in Figures 2A-2D, and in 0.1 mM potassium phosphate at pH 7.2 and 22°C in Figures 2E-2F. 



data. In the first step, Figure one of JMcMinn and Ott- 
away [20] is used to determine the parameters Kf„A> 
KmB and K^c- Because these data describe the 
changes in the initial rate of OGDHC activities in 
response to variations in the concentration of differ- 
ent substrates in the absence of products, they faciU- 
tate identifying the kinetic parameter associated with 
the binding of 2-oxoglutarate, CoASH*', and NAD' . 
Next, the kinetic data from Figure three of McMinn 
and Ottaway [20] are used to estimate Kia and K^r- 
Finally, the kinetic results from Figure five of Smith et 
al. [10] are used to determine the last two parameters 
related to Succinyl-CoA*' inhibition: Ktc and q. 



Henceforth, these kinetic parameters are fixed at their 
estimated values (Table 1). 

Parameterization of the cofactor-dependent regulatory 
mechanisms 

Denton et al. [23] conducted a number of experiments 
to study the effects of Ca^*, pH, and adenine nucleo- 
tides on the activity of OGDHC from pig heart mito- 
chondria. Their data are used here to identify the 
kinetic parameters that characterize the activation/ 
inhibition mechanism of Ca * , pH, and adenine 
nucleotides (See Figure caption for details). Enzyme 
activities are expressed as units of enzyme activity per 
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Figure 3 OGDHC activity regulated by calcium and pH. (A) OGDHC activity as a function 2-oxoglutarate concentration with 1.0 mM NAD, 
0.25 mM CoASH, and in the presence of 5 mM EGTA (circle) or 5 mM EGTA plus 5 mM CaCI2 (up-triangle)Kinetic data were obtained from 
Figure Three of [23]. (B) OGDHC activity as a function free calcium concentration in the presence of 1.0 mM NAD, 0.25 mM CoA and 0.1 mM 2- 
oxoglurate. Kinetic data were obtained from Figure Four of [23]. (C) OGCHC activity as a function of pH in the presence of 1.0 mM NAD, 0.25 
mM CoASH with additions of 0.1 mM 2-oxoglurate plus 5 mM EGTA plus 5 mM CaCl2 (up-triangle) or plus 5 mM EGTA (circle), or with additions 
of 25 mM 2-oxoglutarate plus 5 mM EGTA plus 5 mM CaCl2 (square) or plus 5 mM EGTA (down-triangle). Kinetic data were from Table 2 in [23]. 
Solid lines are model fitting results for the data points represented by symbols. Measurements were conducted at pH 7.0 and 30°C for Figure 3A 
and pH 6.8 and 30°C for Figure 3B and all in 50 mM MOPS buffer plus 1 mM dithiothreitol and with additions of appropriate substrates. 



mg of protein. (One unit of activity is the amount of 
enzyme which transforms 1 |imol of substrate per min- 
ute at 30°C). 

To fit the data that describe the regulatory effects of 
Ca^^ from Figures 3 and 4 of Denton et al. [23], it is 
necessary to account for mechanisms of allosteric 
activation/inhibition of Ca * in the model. It has been 
shown that Ca^^ can significantly affect the function 
of OGDHC by modulating the enzyme affinity for 2- 
oxoglutarate, showing a sigmoidal kinetics. Here we 
propose a general scheme of nonessential and mixed- 
type activation to characterize the effects of Ca 
Based on this scheme, we assume the presence of two 
binding sites for Ca^* on OGDHC and modify Vmax 
and Km of the enzyme complex for 2-oxoglutarate as 
follows: 



= V, 



max,0 



1 + -I- J 



1 + 



[Ca^ 



OlCaKaCa 
^2 4 



[Ca^^]^ 

2+ .2 



/ [Ca^M [Ci^n \ 
1 + + T 



(14) 



[Ca^ + ] [Ca^ + ] 



1 + 



The modified flux expression for OGDHC is obtained 
by substituting Equation (14) into Equation (11). Three 
adjustable parameters {(Xcm Pca> and KaCa) are estimated 
based on the data from Denton et al. [23]; the model 
fits are shown in Figure 3A-B. For model simulations, 
the kinetic constants for substrates: 2-oxoglutarate, 
NAD', and CoASH*' are fixed at their previously 




0.6 0.9 1.2 
[ADP] or [ATP] (mM) 



2 3 
[ADP] or [ATP] (mM) 




1 1.5 2 
[ADP] or [ATP] (mM) 



Figure 4 OGDHC activity regulated by nucleotide with and without calcium (A, B) OGDHC activity as a function of various concentrations 
of ATP (circle), ADP (up-triangle) or ADP with 1.5 mM ATP (down-triangle) in the presence of 1.0 mM NAD, 0.25 mM CoASH, with the addition 
of either 0.1 mM 2-oxoglurate plus 5 mM FGTA plus 5 mM CaCI2 (A) or 2 mM 2-oxoglurate plus 5 mM EGTA (B). Kinetic data were obtained 
from Figure Five of [23]. (C) OGDHC activity as a function of ADP/ATP ratio in the presence of 1.0 mM NAD, 0.25 mM CoASH and either 0.1 mM 
2-oxoglurate plus 5 mM EGTA plus 5 mM CaCl2 (circle) or 2 mM 2-oxoglurate plus 5 mM EGTA (up-triangle). Kinetic data were obtained from 
Figure Six of [23]. Solid lines are model fitting results for the data points represented by symbols. 
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estimated values, obtained from data of McMinn and 
Ottaway [20]. Based on the fits to these data we find 
that Pca is close to one. For simplicity, we fixed Pca = 1. 
meaning that Ca^* affects only the of 2-oxoglurate, 
not the Vmax- This mechanism is also consistent with 
the conclusions of Denton et al. [23]. (Statistical analysis 
of different model formulizations supports the validity 
of the null hypothesis that adding the extra adjustable 
parameters, Pcw does not lead to any significant 
improvement in fitting results.) The estimated values of 
the kinetic parameters are summarized in Table 1. 

Experimentally it has been shown that the maximal 
activity of OGDHC is largely unaffected by changes in pH 
over the range 6.6-7.4, whereas the of the enzyme is 
markedly altered by pH in this range [23]. In our model, 
the effect of pH on the OGDHC activity was described 
based on the observations studies of Denton and collea- 
gues [23] . Here, protons are treated as the essential activa- 
tors of OGDHC which increase the binding affinity of the 
enzyme to 2-oxoglutarate. Therefore, the of 2-oxoglu- 
rate is modified by multiplying the term /^///[H^] such 
that I<^A_ 1 = {K^A. i^<a//)/[H*] (in Equation (15)). Figure 
3C illustrates the model fits to the data obtained from 
Table 2 in McCormack and Denton [23] where the activity 
of ODGHC was studied under varying pH in both pre- 
sence and absence of Ca ^ in the buffer. 

Denton and colleagues [23] also studied the effect of 
adenine nucleotides (ATP and ADP) on the OGDHC 
activity where it was shown that both ATP and ADP 
significantly impact the of the en-zyme for 2-oxoglu- 
tarate and that the regulations of OGDHC by Ca^^ and 
adenine nucleotides seem to be independent. Here the 
regulatory effects of nucleotides on OGDHC activity is 
modeled as similar to that of Ca^*. Specifically, we 
assume that there are different binding sites on the 



OGDHC that bind to ATP or ADP. (The available 
kinetic data cannot exclude the other possibility that 
ATP and ADP can bind at the same site.) Therefore, the 
Vjnax and for 2-oxoglutarate was modified as a func- 
tion of nucleotide concentrations as follows: 



V oiAwKAri' J V 



,.|ADP|i 



I ^ [ATPIt \ / ^ |ADP|t \ 
[ATPIt 



CADI'KiADI' , 

[ADPIt , 



KmA,2 = K,„A,l 



/ ^ [ATPItA / ^ [ADPjrX 

\ KjAW ) \ KlALW / ^<lH 

/ ^ JArP|^\ / ^ JADP]^\ |H-] 



(15) 



where [ATP]t and [ADP]t represent the total concen- 
trations of the nucleotides. Experimental results show 
that neither free nucleotides nor magnesium binding 
nucleotides are solely responsible for activation 
observed. More data are needed to quantitatively and 
qualitatively specify the activation effects of each nucleo- 
tide's ionic forms. Figure 4(A-C) are model fits to the 
data obtained from Figures 5 and 6 of Denton and col- 
leagues [23]. Optimization results based on these data 
indicate that both ATP and ADP change the K,„ of 
OGDHC for 2-oxoglutarate without altering the maxi- 
mum activity, which suggests that both Patp and Padp 
are equal to one. 

Mg^^ is known to regulate the activity of OGDHC. In 
a recent study, Rodriguez-Zavala et al. [37] examined 
the effects of ligands, such as ATP, ADP, Ca^"^, and Mg^ 
* on the activity of OGDHC in both isolated pig heart 
enzyme complex and mitochondrial extracts. These data 
facilitate the characterization of the regulatory effect of 
Mg * on the OGDHC activity and are used here to esti- 
mate the Mg * associated kinetic parameters. Enzyme 




[ADP] (nM) [aKG] (mM) [aKG] (mM) 



Figure 5 OGDHC activity regulated by magnesium (A) OGDHC activity of isolated enzyme as a function of ADP in the presence of 1.0 mM 
NAD, 0.25 mlVl CoASH, and 0.5 mM 2-oxoglurate with additions of 0.0 (circle), 25 (jM (up-triangle), 50 |jM (down-triangle) and 200 |jM (square) 
Mg^*. Kinetic data were obtained from Figure Two of [37]. (B-C) OGDHC activity as a function of 2-oxoglutarate in the presence of 1.0 mM NAD, 
0.25 mM CoASH with 600 [iM (B) or 0 pM (C) Mg^* addition, (up-triangle) no further additions (control); (circle) 4-600 pM ATP; (down-triangle) 
4-500 pM ADP. Kinetic data were obtained from Figure Four of [37]. V^nax has been changed less than 10% to better fit the symbols. Solid lines 
are model fitting results for the data points represented by sym-bols. The experimental data were obtained at pH 7.35 and 30°C in 120 mM KCI, 
20 mM MOPS-K 0.5 mM EGTA buffer 
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log[Ca*'l (M) [ADP]/[ATP] log([NADH]/[NAD]) 



Figure 6 Predictions of OGDHC activity comparing with experimental results. (A) OGDHC activity as a function of free calcium 
concentration in present of 2.0 mM NAD and 0.25 mM CoASH in the presence of either 1.5 mlVl ATP (circle) or 1.5 mM ADP (up-triangle). 
Concentrations of 2-oxoglutarate were 0.05 mM (ADP alone) and 0.1 mM (ATP alone). Kinetic data were obtained from Figure Three of [38]. (B) 
OGDHC activity as a function of ADP/ATP ratio in the presence of 2.0 mM NAD, 0.25 mM CoASH and 0.05 mM 2-oxoglutarate in the presence of 
either 1.0 nM (circle) or 1.0 |jM (up-triangle) free Ca^*. Kinetic data were from Figure Four of [38]. (C) OGDHC activity as a function of NADH/ 
NAD ratio 0.35 mM Mg'^* and 0.4 mM total nicotinamide nucleotide in the assay. Ca"^* concentration was either 1 nM (circle) or 0.1 mM (up- 
triangle) and 2-oxoglutarate concentration was 10 mM (circle) or 0.375 mM (up-triangle). Kinetic data were obtained from Figure Five b of [38]. 
Solid lines are model predictions for the data points represented by symbols. The experimental data were obtained at pH 7.2 and 30°C in 
reaction medium of 50 mM MOPS, 75 mM KCI, 60 mM sucrose, 2 mM KH2PO4, 1 mM EGTA. 



activity is measured in nmol NADH^^ formed per min- 
ute per mg protein. 

Experimental data from purified OGDHC from pig 
heart mitochondria from Rodriguez-Zavala et al. [37] 
with zero Mg^* were used to estimate the V^^^^ of the 
enzyme. Figure 5B shows the model simulations (lines) 
using the parameter estimates obtained above, with the 
exception that Vmax was adjusted to match these experi- 
mental data. Data in Figure 5A with non-zero Mg^^ con- 
centration and the data in Figure 5C are then used to 
estimate the kinetic parameters associated with the bind- 
ing of Mg^*. Our model fits to these data assume that 
Mg * not only increases the activity by binding to the 
enzyme complex, but also potentiates the 2-oxoglutarate 
affinity to the enzyme and decrease the K,^ of OGDHC 
for 2-oxoglutarate. Specifically, we assume two binding 
sites for Mg^* and modify the Vmax and the Km as 

,2 + ,2 



V, 



max, 3 ~ *'max,2 



^ , [Mg^M , [Mg- 



2+l2 



V aMgK.Mg {aMgKaMg) ) 

i , [Mg2^] [Mg2-f \ 
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(16) 
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The data shown in Figure 5 are used to estimate the 
adjustable kinetic parameters related to Mg^^ ions in 
our kinetic model for OGDHC. Fits to the data are 
plotted in Figure 5 and the parameter values 



summarized in Table 1. The developed model is able to 
satisfactorily explain the effect of Mg * ions on the 
enzyme activity. These results, combined with those 
shown in Figures 2 and 3, imply that the matrix free 
Ca^* and Mg^* ions concentrations exert significant and 
distinct effects on the OGDHC activity. 

Complete flux expression for the 2-oxoglutarate 
dehydrogenase complex 

Based on the proposed mechanisms of allosteric activa- 
tion and inhibition of various cofactors, the flux expres- 
sion (Equation (11)) of the OGDHC can be further 
modified. Applying the catalytic and regulatory mechan- 
isms of Equations 14-16, the final forward flux expres- 
sion is 



V/[A][B][Cj-N 



ogdiic 



lAl[Bj[C] + ^i^„c[A][B][R] 

Kir 

+J^mclAl[Bj + -^K^b[A][C][Q] 
+J^™b[A][C] + :^^^[A][Q][R] 



KiqKir 



+ — i<^mc|Bj[R] + [QIIR] 



\ Kir 



KiqKr 



(17) 



[Ca2 * I |Q 



where 



lATP].,. 



|H*| 



Qi et al. BMC Biochemistry 201 1, 12:53 
http://www.biomedcentral.eom/1 471-2091/1 2/53 



Page 11 of 1 5 



\ (^MgKaMg {aMgKaMgY/ 

Estimated values of KaH> Kiatp, Kuadp, KaMg> cCadp> 
CCatP' Olca> OlMgi and f^Mg are listed in Table 1. 

Independent validation of the developed kinetic model of 
OGDHC 

Finally, the model is independently validated (corrobo- 
rated) by comparing the model predictions to the initial 
rate data of Rutter and Denton [38] on the kinetics of 
OGDHC obtained from permeabilized mitochondria and 
mitochondrial extracts (see Figure 6(A-C)). They studied 
the regulations of NAD-linked isocitrate dehydrogenase 
and 2-oxoglutarate dehydrogenase by Ca nucleotide 
and nicotinamide nucleotides in permeabilized rat heart 
mitochondria and in mitochondria extracts. Data from 
their study were not used for estimation of model para- 
meters and used here to further validate the proposed 
mechanisms and regulation of OGDHC. Therefore, the 
flux expression of Equation (17) was used for simula- 
tions with the values of the kinetic parameters the same 
as estimated before (see Table 1). The model accurately 
describes the kinetics and regulation of OGDHC, 
observed experimentally, without having to re-estimate 
the model kinetic parameters, signifying the accuracy of 
the model and the associated model parameters. 

To determine the degree to which the model simula- 
tions are sensitive to the estimated parameter values, 
the relative sensitivities are computed and listed in 
Table 1. A high sensitivity value indicates that a small 
change in a given parameter can lead to significant 
changes in model outputs, used to identify the para- 
meter values. All of our adjustable parameters of the 
model have sensitivities over 30%. Two parameter esti- 
mates (KaCa and KaMg) show relatively low sensitivity 
compared to the others, indicating that predictions of 
the developed model are less sensitive to these two 
values. This implies that these two parameters may not 
be identified accurately by the present analysis, given 
the sparseness of the data sets analyzed in this work 
regarding the regulation of OGDHC by Ca^^ and Mg^^. 
Further experiments are required to adequately establish 
the appropriate regulatory mechanisms and the robust- 
ness of each model parameters. 

Discussion 

A number of kinetic models have been previously devel- 
oped to explain the basic catalytic mechanisms and reg- 
ulations by cofactors of OGDHC. Sanadi et al. [12] first 
proposed a Hexa-Uni-Ping-Pong mechanism for the 



overall reaction by studying various roles and locations 
of the cofactors: thiamine pyrophosphate, lipoic acid, 
and FAD^' within this complex. Hemada et al. [19] con- 
ducted kinetic studies and proposed a similar mechan- 
ism to that of Sandi et al. and suggested that NADH 
is a competitive inhibitor of NAD'. Whereas Smith et 
al. [10] suggested a noncompetitive inhibition of 
NADH ' with NAD , the catalytic mechanism was not 
consistent with Sandi et al. Later, McMinn and Ottaway 
[20] tested a series of possible alternate mechanisms 
using computer optimization techniques and initial velo- 
city studies and concluded that the binding of NAD' 
and CoASH ' and the release of Succinyl-CoA ' is a 
random order, whereas the binding of the substrate 2- 
oxoglutarate and release of the product CO2 still follows 
a Ping-Pong mechanism. Besides above experimental 
studies, a number of integrated models of mitochondrial 
bioenergetics have been developed which used different 
type of OGDHC models. Cortassa et al. [39] describe 
the activity of OGDHC as a function of Ca^"^, Mg^"^ and 
substrate concentrations using phenomenological terms. 
Wu et al. [40] used a simple kinetic model of OGDHC 
from Kohn and Garfinkel [41] in their integrated model 
of TCA cycle that does not incorporate the regulatory 
effect of metal ion cofactors. In a recent integrated 
study of mitochondrial bioenergetics, Bazil et al. [42] 
developed a kinetic model of OGDHC based on a Hexa- 
Uni-Ping-Pong mechanism with a general description of 
the cofactor dependency of OGDHC activity. In sum- 
mary, there have been a wide variety of kinetic models 
of OGDHC with contrasting kinetic mechanism and 
cofactor regulations. 

In this paper, we developed a unified mechanistic 
model of OGDHC, in which Ca^+, Mg^^, ADP^', and pH 
are treated as activators and ATP*' as inhibitor of the 
OGDHC activity. The present model offers more realis- 
tic and meaningful explanations on the catalytic proper- 
ties and regulation mechanisms of OGDHC than 
previous attempts. The analysis also provides a unique 
set of kinetic parameters that consistently describe a 
wide variety of experimental data sets on OGDHC func- 
tion, obtained from diverse sources. Based on the 
assumed ter-ter mechanism and associated conforma- 
tional changes, we are able to consistently reproduce the 
observed kinetics of OGDHC with a minimal number of 
model parameters. Thus, the proposed mechanism is 
found to be more appropriate compared to other alter- 
nate kinetic models [43]. 

Effects of nucleotides on the OGDHC activity 

Energy-linked regulators, ADP and ATP, as well as inor- 
ganic phosphate, have been investigated for over two 
decades for their profound effects on kinetic properties 
of OGDHC. Kinetic studies of mammalian OGDHC, 
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isolated from varied sources, have shown that ADP 
causes activation of OGDHC [44-47]. This enzyme com- 
plex is sensitive to ADP, where ADP significantly 
decreases the for 2-oxogluterate without affecting 
the maximum rate (V^^^) of the reaction via allosteric 
interactions. For example, studies on the OGDHC of rat 
heart mitochondria show a seven fold decrease in the 
Km value for 2-oxogluterate by ADP and thereby 
strongly increases the affinity of OGDHC for the sub- 
strate [47]. Other studies of OGDHC from human heart 
make similar conclusions on the activating effect of 
ADP [48]. It has also been shown that, at subsaturating 
concentration of 2-oxoglutarate the relationship between 
initial reaction rates of OGDHC and concentration of 
ADP is sigmoidal, suggesting a positive cooperativity in 
binding of ADP to the enzyme complex [49]. In con- 
trast, at a suboptimal 2-oxoglutarate concentration, ATP 
was shown to inhibit OGDHC activity in pig heart and 
bovine kidney mitochondria [23,45]. In addition direct 
inhibition effect, recent investigations have shown the 
possible indirect inhibition OGDHC by ATP because of 
chelation of divalent ions which activate OGDHC, such 
as Ca^"^, Mg^"^ [37]. The activating action of ADP and 
inhibiting action of ATP are in competitive opposition 
(Figure 3E). It is still unclear if these two effectors bind 
on the same site on the complex or not. Model analysis 
based on the available kinetic data cannot exclude either 
possibility. In our mechanistic model, the regulatory 
effect of ADP and ATP is incorporated by assuming dif- 
ferent binding sites for ADP and ATP in the enzyme 
complex and the model satisfactorily describes the acti- 
vating effect of ADP and inhibitory effect ATP observed 
in many experiments [46]. 

While Zavala et al. interpret the data of Figure 5A to 
indicate that MgADP' is the effective activator of 
OGDHC activity, our model analysis of the available 
data sets on Mg ADP, and ATP dependent kinetics 
(Figures 4 and Figure 5B and 5C) reveals that, magne- 
sium and ADP have independent parallel effects on the 
OGDHC activity, the most parsimonious explanation of 
the data. However, al-ternative, more complex, models 
cannot be ruled out. 

Like ADP, Pi has also been shown to decrease the 
value for 2-oxogutarate, without affecting the V^^^ of 
OGDHC reaction [44,50]. In a recent report, the Pi acti- 
vation showed biphasic behavior, with pH dependence 
[37]. In the physiological concentrations range Pi exerts 
monophasic activation of OGDHC [37], which can be 
descried by Equation (15) with three extra parameters. 
Due to the lack of consistent kinetic data for Pi effects, 
we do not integrate a Pi dependent regulation mechan- 
ism in our current model. However, this energy linked 
effector may be physiologically important. The overall 
rate of oxidative phosphorylation is largely determined 



by phosphorylation potential [51,52]. In cells when ATP 
utilization increases, the production of ADP and Pi 
increase. Therefore, activation of OGDHC by ADP and 
Pi may represent a compensating effect. 

Effect of pH on the OGDHC activity 

Mitochondrial matrix proton (H* ion) concentration is 
known to affect the OGDHC properties. Specifically, 
studies on pig heart OGDHC showed that the change in 
pH in the range 6.6 to 7.4 can significantly alter of 
the enzyme for 2-oxogluterate, without affecting its 
maximal activity. McCormack and Denton illustrated 
the effect of pH on OGDHC activity both in the absence 
and presence of Ca^^ in their assay mediums [46]. Other 
experimental observations have shown that hydrogen 
ions favor the higher affinity of OGDHC for 2-oxogluta- 
rate [45,53]. In our model we hypothesize that hydrogen 
ions are essential activators of OGDHC activity to 
describe the observed pH dependency of the OGDHC 
kinetics. 

Effects of Ca^\ Mg^^ and EGTA on the OGDHC activity 

Studies by McCormack and Denton demonstrate the 
activating effects of Ca^* ions on intra-mitochondrial 
dehydrogenases: pyruvate (PDH), NAD-isocitrate (NAD- 
ICDH), and 2-oxoglutarate (OGDHC) [24,54,55]. Speci- 
fically, the rise in cytosolic Ca^* concentration in 
response to extrinsic stimuli, such as hormones can 
enhance mitochondrial oxidative metabolism via direct 
activation of these three Ca^* sensitive dehydrogenases. 
Such mechanisms may serve as a complementary way to 
stimulate ATP-synthesis to meet the increased energy 
demand of the cell [24,54,55]. 

Mg^^ ion has also been shown to regulate the 
OGDHC activity either by directly activating the enzyme 
or by modulating the Ca^^ effects on the enzyme. In 
some studies, Mg ^ shows no effects on OGDHC activ- 
ity [23]. However, in other studies, Mg ^ is shown to 
increase the maximal activity of the enzyme and the affi- 
nity of OGDHC for 2-oxogluterate by enhancing the 
Ca^* stimulatory effects on the enzyme complex 
[9,25,37]. These different observations could be 
accounted for the different levels of endogenous Ca^* 
and Mg * present in the purified enzyme complex pre- 
pared by different methods. Another possible explana- 
tion is that the stimulatory effects of Mg^^ is TPP- 
dependent, which is not explicitly considered in our 
model. Panov and Scarpa [9] found that Mg * only 
exerts its stimulatory effects in the presence of TPP, 
though exclusion of TPP from the reaction medium has 
no effect on the initial enzyme activity in the absence of 
Mg^^. Also, it has been clearly shown that Mg^^ may 
affect the rate of oxidative phosphorylation in isolated 
mitochondria primarily via modulating the OGDHC 
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activity [25]. The site of action of Mg^+ ion on OGDHC 
is unknown. In the present model, we hypothesized a 
general scheme of nonessential activation of Ca^*, by 
considering two Ca^^ and Mg^* binding sites on 
OGDHC. The Mg ^ effect is incorporated in our model 
by exclusively modifying the enzyme activity and 2-oxo- 
glutarate binding step. So the parameters Knax and 
are accordingly expressed as functions of Mg * (Equa- 
tion (16)). Currently, our model assumes that the turn- 
over rate of El is modified to same value for binding 
either one ion or two ions. And to make it simple, our 
model does not include possible interaction between 
Mg * and Ca * at high concentrations either. The effects 
are not additive [9] at high concentration, suggesting 
that Mg^^ and Ca^^ may compete for the binding site. 
Additional kinetic data set are necessary to test different 
mechanisms and refine our model to more accurately 
describe the nature of cation dependent kinetic of 
OGDHC. 

EGTA, which is used in many studies to control Ca^* 
ion concentration in reaction media, has been shown in 
experiments and theoretical analysis to inhibit the 
NAD-linked isocitrate dehydrogenases (ICDH) through 
the binding complex, MgEGTA [36]. To date, it is still 
not clear if there is similar inhibition effect of EGTA or 
EDTA on the activity of OGDHC. McCormack and 
Denton [23] concluded that the OGDHC sensitivity 
toEGTA is very similar to that observed with ICDH 
[56], because addition of calcium chelators EGTA or 
EDTA is associated with a marked decrease in the activ- 
ity of OGDHC at 0.2 mM 2-oxoglutarate. Panov and 
Scarpa [9], in ascribing the inhibition effect of EGTA to 
the complex formation between Ca * and chelators, 
concluded that the effected of Ca * and chelators is 
associated with different endogenous cation levels in dif- 
ferent preparations. But this explanation cannot account 
for McCormack and Denton's observation that EGTA 
or EDTA causes a 40% decrease of activity of OGDHC 
after using Chelex remove much of the endogenous Ca 
^ in the buffer. It is also noted that the Km for 2-oxoglu- 
tarate in the absence of Ca^* is 4 ± 1.1 mM measured 
by Panov and Scarpa [9] for commercially available 
enzyme (Sigma, St. Louis, lot 44H80801), which is 
almost 15 times the estimated value used in our model 
(Table 1) and that found by previous workers [19,27]. 
Only by using the reported Km of Panov and Scarpa [9], 
can we reproduce their data using the same mechanism 
(Equation (17)). In the absence of clear experimental 
evidence and sufficient data set, our model does not 
explicitly account for an inhibition effect of EGTA or 
EDTA. 



Conclusion 

Our mechanistic OGDHC model based on a detailed 
catalytic mechanism successfully provides a single con- 
sistent theoretical explanation for many previously unre- 
solved experimental observations on the kinetics and 
regulations of OGDHC. In particular, it suggests the 
most plausible physiologically regulations of OGDHC by 
NAD(H) oxidation-reduction state, the nucleotide phos- 
phorylation potential, pH and various metal ions (Mg^^ 
and Ca^*). As a rise in NADH can reduce the OGDHC 
flux and thereby provides feedback regulation through 
the electron transport chain, it is important to ask how 
NAD oxidation-reduction state and oxidative phosphor- 
ylation state exert a coherent regulation of OGDHC in 
physiological context. Furthermore, how does the 
OGDHC respond to stimuli via the mitochondrial Ca ^ 
transport system? Such questions may be addressed by 
applying the present model in an integrated framework 
[43] along with other dehydrogenases [36], the oxidative 
phosphorylation system [57], electron transfer system 
[58], and cation transport systems [59-61]. 
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